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© ProcedS de demodulation coherente d'un signal module numeriquement en phase continue et a 
enveloppe constante. 

© Le procedS assure la demodulation coherente par traitement numerique d'un signal module en phase 
continue (par exemple de type OMSK). Le signal regu est transpose en bande de base, convert en numerique 
<et transmis au processeur de signal. Chaque paquet de reformation binaire transmise comporte une sequence 
_preliminaire connue sur N bits permettant I'estimation approximative de synchro-trame et synchro-bit d une part, 
Sde la phase initiale et de I'ecart residue! en frequence d'autre part. L'affinement progressif desdites estmatons 
©est obtenu au moyen de deux boucies numeriques imbriquees : une boucle lente pour la detection de synchro- 
— bit et une boucle rapide effectuant des decisions intermediates sur des blocs additionneis de bits pour 
^'estimation de la phase initiale et de l'£cart residuel en frequence. 
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ProcedS de demodulation cohtSrente d'un signal moduli numeriquement en phase continue et a 

enveloppe constants 

L'invention conceme un proc§de de demodulation coherente par traitement numerique d'un signal 
module numeriquement en phase continue et a enveloppe constante, le terme module de ladite phase etant 
ggal au produit de convolution de rimpulsion de phase etalee sur plusieurs temps bits par rinformation 
binaire transmise par paquets, le signal regu transpose en bande de base sur deux voies en quadrature 
5 etant convert! en numerique et transmis a un processeur de signal qui effectue le traitement du processus 
de demodulation. 

Ce proc§d«§ peut s'appliquer a n'importe quelle modulation de phase du genre precite : (GMSK, MSK, 
2SRC, TFM, GTFM...) dont la loi devolution suivant une variation progressive de la phase a pour avantage 
de reduire le spectre. Un autre avantage r§sulte du fait que l'§nergie transmise est constante. 
10 En particulier la modulation de type GMSK dont la variation de phase est etalee sur (e plus grand 
intervalle de temps (5 temps bit), presente la meiileure efficacite spectraie. Malheureusement, ceci a pour 
effet d'augmenter d'une facon importante Interference inter-symbole. 

Uutilisation de cette modulation a bande etroite peut s'envisager dans de nombreux domaines tels que 
les systfcmes de communications protegees en VHF et UHF, les transmissions par satellites ou les reseaux 
rs radiomobiles. Les avantages precitSs I'ont fait retenir par le Groupe Special Mobile (GSM) du CEPT pour 
§tre utilisee dans le futur reseau mobile numerique Pan-Europeen a partir de 1992. 

Les divers procedes de demodulation connus mettent en oeuvre des methodes differentieile ou 
coherente. 

La premiere methode a Tavantage d'etre relativement simple mais les performances en terme de taux 
20 d'erreur sont tres d^gradees. 

La demodulation coherente presente de meilleures performances mais elle necessite un dispositif 

supplemental de recuperation de la phase de la porteuse. 

Une des faiblesses rencontrees avec ce type de demodulation reside dans Pemploi de methodes 
classiques de synchronisation qui utilisent des boucles a verrouillage de phase pour recuperer la porteuse 
25 et Phorloge. 

En effet dans le cas d'un systeme fonctionnant en Acces Multiples a Repartition dans le Temps 
(AMRT) ou en evasion de frequence (EVF) et lorsque le signal subit des evanouissements dus au canal, les 
temps de ^synchronisation des boucles analogiques deviennent trop longs et reduisent la duree utile du 
signal (cf. brevet americain N* 4 570 125 de R.B.Gibson et B.Hill). 
30 La mise en oeuvre d'une methode de demodulation coherente par traitement numerique du signal offre 
comme principal avantage la possibility de memoriser et de traitor le signal par paquets pour chacun 
desquels on doit effectuer un traitement sequential se terminant par une decision sur les informations 
binaires transmises. 

La premiere operation de la sequence consiste a rechercher le debut du paquet ; c'est la synchro- 
35 trame. Puis la synchro-bit permet de determiner les instants de decision et d'assurer le bon caiage 
temporel du filtre adapts. Le role de ce filtre est de reduire le bruit sans degrader ['information utile. 

Le dernier traitement est extremement important : il s'agit de ('estimation de la phase initiate et de 
l'£cart residuel de frequence. 

La phase initiale est un parametre qui n'est pas maftrise dans un systeme de transmission. Une 
40 mauvaise estimation de ce parametre est desastreuse sur le taux d'erreurs. 

L'Scart de frequence residue! est la resultante de I'ecart de frequence entre I'emetteur et le recepteur et 
de Pecart en frequence dQ a I'effet Doppler. Une evaluation defectueuse de cet ecart de frequence se 
traduit par des erreurs de decision sur les derniers bits du paquet, lorsque la phase a suffisamment tourne 
pour provoquer de telles erreurs. 
45 Aprfcs avoir estime la phase initiale et Tecart de frequence residuel on effectue une compensation et 
finalement on decide les bits transmis. 

Une mdthode numerique de demodulation utilisee en 2SRC a ete proposee dans Particle de LOUBA- 
TON et VALLET intitule : " Demodulation pseudo-coherente de signaux de type MSK adaptee aux 
transmissions en EVF" et publie dans la Revue Technique Thomson-CSF, vol. 17, Septembre 1985, N 3, 
so pages 521 -554. 

Dans cette methode on retrouve la sequence de traitement suivante : synchro-trame par correlation 
partielle ; synchro-bit par detection de passage a zero de la phase drfferentielie ; filtrage adapte ; estimation 
de Tecart de frequence residuel par transformation de Fourier rapide sur les echantillons au carre ; 
estimation de la phase initiale par moyenne ; compensation en phase. 

2 
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Apres simulation, il s'avSre que les algorithmes proposes s'adaptent mai k la GMSK. 
En effet. k cause de I'interference inter-symbole qui ne peut plus etre negligee, la determination de la 
synchro-bit est trfes degradee en presence d'un €cart residuel en frequence superieur k 200 Hz. 

De plus pour estimer cet ecart de frequence, on ne peut plus eliminer la modulation en effectuant une 
5 elevation au carre. 

Une technique de demodulation pour la transmission de paquets par radio est donnee dans Particle de 
C.HEEGARD, J.A. HELLER et AJ. VITERBI intitulS : "A microprocessor-based PSK Modem for Packet 
Transmission over Satellite Channels" et paru dans IEEE, vol.COM-26, N* 5, Mai 1978, pages 552 a 564. 
Inspire de cette technique, qui ne s'applique qu'aux modulations de type PSK sans interference 

10 intersymbole, le procede de ^invention vise a obtenir une synchronisation permettant d'effectuer la 
demodulation coherente de n'importe quel type de modulation presentant de IMnterference intersymbole et 
ce, meme avec un niveau de bruit et un ecart residuel en frequence importants. 

A cet effet ce procede est remarquable en ce que chacun desdits paquets comporte une sequence 
preliminaire par laquelle on dispose d'un signal de reference connu sur N bits permettant de detecter 

75 approximativement la synchro-trame et la synchro-bit par correlation avec la phase differentielle et 
d'amorcer I'estimation des parametres de phase initiale Bo et d'ecart residuel en frequence Afo, I'affinement 
progressif des valeurs approchees etant obtenu au moyen de deux boucies numeriques imbriquees : une 
boucle lente initialisee au-dela d'une valeur de seuil pour la detection de la synchro-bit, une boucle rapide 
effectuant des decisions intermediaires sur des blocs additionnels de bits pour i'estimation de do et Afo. 

20 Lesdites detections de synchro-trame et de synchro-bit sont obtenues par une premiere correlation sur 
la phase differentielle permettant de connattre I'instant d'emission desdits paquets a ± T/4, T etant la dures? 
d'un bit, puis par une deuxieme correlation sur la phase differentielle effectuee avec un signal de reference 
decaie de T/4. Les fonctions de correlation correspondantes presentent chacune un pic independant de la 
phase initiale et tr§s peu dependant de I'ecart residuel en frequence, le pic de niveau superieur et le pic de 

25 niveau inferieur definissant respectivement une synchro-bit principale SYNP et une synchro-bit secondaire 
SYNS. La precision de ± T/8 ainsi obtenue sur la synchro-bit est suffisante pour connattre i'instant 
d'echantillonnage. 

Ladite detection de synchro-bit est suivie d'un fiitrage adapte effectue avec la valeur SYNP par un filtre 
a reponse impulsionnelie finse de type Gaussien afin de limiter la bande de bruit 
30 Lesdites estimations de la phase initiale So et de I'ecart residuel en frequence Af 0 a la suite dudit 
fiitrage comportent les etapes suivantes : 

- Elimination du terme de modulation en effectuant le produit du signal regu par le conjugue du signal de 
reference. 

- Derouiement de la phase en eliminant les sauts de phase de 2w pour obtenir une variation lineaire ayant 
35 pour equation y = A«ox + do avec A«o = 2-nAfo. 

- Calcul des parametres estimes Awo et ?o par une methode de regression lineaire et de I'ecart 



= Wr <Awo.X,+ Go) ' 



40 i=0 

entre les points correspondant a ladite phase deroulee et ladite droite de regression. 

Si ledit ecart e est inferieur a ladite valeur de seuil, ladite estimation de Afo et Qq est affinee suivant une 
boucle rapide en plusieurs passes exploitant les decisions intermediaires sur les N bits de la sequence 
preliminaire auxquels on ajoute a chaque passe un certain nombre de bits decides. 

Si ledit ecart e depasse ladite valeur de seuil du fait d'une evaluation defectueuse de synchro-bit. ledit 
processus de calcul est reinitialise suivant une boucle lente pour refaire le fiitrage adapte et I'estimation de 
Afo et do k partir de i'autre valeur de synchro-bit egale k ladite valeur secondaire SYNS. 

Apres la derniere passe, une compensation ne laisse subsister que la composante de phase du signal 
qui n'est plus affectee par I'ecart residuel en frequence ni par la phase a I'origine. 

La decision finale est ensuite effectuee puis un decodage differentiel fournit enfin la suite d'informations 
binaires transmises. 

L'invention sera mieux comprise a I'aide de la description suivante donnee a titre d'exempie non 
limitatif, ladite description etant accompagnee de dessins qui represented : 
Figure 1 : le schema bloc d'un dispositif modulateur-demodulateur. 

Figure 2 : les variations d'impulsion de phase pour les modulations de types GMSK, MSK et 2SRC. 
Figure 3 : les diagrammes d'occupation spectrale pour les modulations de types GMSK, MSK et 
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2SRC. 

Figure 4 : le diagramme de I'oeil pour !a modulation GMSK. 

Figure 5 : les variations temporelles de la phase du signal regu apres filtrage lors de la sequence 
devaluation de Oo et Afo. 

Rgure 6 : I'organigramme de I'ensemble du traitement de demodulation selon le procede de 

invention. _,.. ei/ 
Rgure 7 : les courbes de taux d'erreur relevees dans la literature pour la modulation GMbK. 
Figures 8 et 9 : les courbes de taux d'erreur selon le procede de demodulation de Invention 
appiiqu£ k un signal moduli en GMSK. 

La demodulation d'un signal GMSK selon le procedS de Invention a ete simuiee sur un dispositif 
modulateur-demoduiateur dont la figure 1 "donne le schema sous la forme de blocs fonctionnels comportant 
successivement : 

- Un ensemble de generation de trame 1 contenant un polynome generateur d'un train binaire pseudo- 
alcSatoire au d<§brt de 16 kbrts/s. Le format de chaque trame transmise est de 128 bits avec une sequence 
preliminaire connue de N = 16 ou 32 bits que Ton vient placer en debut de trame au moyen d'un systeme 
de registres et de bascules. II reste done 128-N bits disponibles pour Information k transmettre. 

- Un modulateur 2 qui g6nere une impulsion de phase a variation progressive de type GMSK. Le signal 
module est disponible en bande de base a partir de deux voies I et Q en quadrature. 

- Les elements de transposition k la frequence intermediate de 70 MHz. Cette transposition est operee a 
remission au moyen des melangeurs 3 et 4 effectuant respectivement le melange des signaux provenant 
des voies I et Q avec le signal d'un oscillateur local 5 a la frequence F e et le meme si gnal dephase de 90 
dans le ddphaseur 6. Apres sommation des signaux en provenance des deux voies dans I'additionneur 7, le 
signal resultant issu de cet addrtionneur traverse successivement un attenuates 8, un generateur de bruit 
blanc gaussien 9 de densiftS spectrale No permettant de sirnuler les conditions reelles de fonctionnement et 

25 un fiitre k large bande 10 centre sur 70 MHz. 

A la reception, le signal transmis est retranspose en bande de base (parties reelle et imaginaire sur les 
voies l' et Q' respectivement) au moyen des melangeurs 3 et 4 , de I'osciilateur local 5 a la frequence F r 
et du d§phaseur 6 . * » . 

- Un ensemble de conversion numerique 11 comportant respectivement pour les deux voies I et Q a 
traiter les filtres passe-bas 12 et 13 qui assurent I'echantilionnage en respectant la condition de Shannon et 
les convertisseurs analogique-numerique 14 et 15 precedes d f echantillonneurs-bloqueurs qui maintiennent 
le niveau du signal pendant la duree de la conversion. 

Les voies I et Q' sont ressorties pour effectuer des controles (entre autres la visualisation des yoies I 
et R apres filtrage) apres transformations inverses operees a travers les convertisseurs numerique- 
analogique 16 et 17 et les filtres 18 et 19. Le train binaire est egaiement sorti apres decodage a travers la 
bascule 20. 

- Un ensemble de traitement 21 comportant un processeur de si gnal dans lequel est effectuee la 
demodulation du signal GMSK selon le procede de ^invention, ce processeur operant en mode complexe et 
etant control^ par un microprocesseur. 

Lorsque Information numerique a transmettre est portee par la phase, le signal module peut se mettre 

sous la forme : 



S(t,B) = V2E exp{j[2^rfot + e 0 + *<t.B)]} 
t : temps 

45 B : (B,) suite d'informations binaires transmises. 
E : energie du signal 
T : duree d'un bit 

fo : frequence de la porteuse (pulsation «o = 27rf 0 ) 
So : phase k I'origine 
so <*>(t,B) : phase variant suivant la suite d' informations binaires 



tp(t,B) = 21f i *Z B. .q(t-iT) 



55 



ou q(t) est I'impulsion de phase de dur6e finie. 

Le terme 1/2 dans I'expression de la phase correspond k IMndice de modulation, e'est-a-dire au rapport 
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excursion de frequence sur frequence rythme. 

La fonction q(t) appelee impulsion de phase traduit la fagon dont la phase va varter. 

La figure 2 represente !a variation de cette fonction q(t) pour les modulations de type GMSK, MSK et 
2SRC. 

5 Pour la GMSK la variation de phase est etalee sur 5 temps bit contre 2 temps bit pour la 2SRC et 1 
temps bit pour ia MSK. 

Cette variation £tant plus lente pour ia GMSK, le spectre occupe est moindre ainsi que le montrent les 
courbes de ia figure 3 representant les variations de la densite spectrale de puissance (DSP) en dB en 
fonction du produit (fT) de la frequence f par la duree T d'un bit, pour les modulations MSK (en trait 
10 continu), 2SRC (en traits pointilles) et GMSK (en traits mixtes). 

Mais I'etalement de {'information sur 5 temps bit pour la GMSK se traduit par la presence d'interference 
inter-symbole mise en evidence par le diagramme de I'oeil de la figure 4 obtenu par I'observation du signal 
GMSK sur un oscilloscope synchronise par Phorloge temps-bit 

On va maintenant expliciter les etapes successives de la demodulation coherente d'un signal module 
15 numeriquement en phase continue et k enveloppe constante selon le procede de {'invention. 

Afin de pouvoir utiliser une methode qui exploite la presence d'interference inter-symbole, ii est apparu 
indispensable d'inserer en tete de chaque paquet une sequence preiiminaire de longueur N. 

Cette sequence connue va permettre de detecter le debut de la trame par correlation et ensuite 
d'amorcer I'estimation de Go et Afo. 
20 Le traitement peut etre decompose en quatre parties principales : synchro-trame et synchro-bit, filtrage 
adapte, estimation de do et Afo, et decision. 



Synchro-trame et synchro-bit 

25 

La detection de la synchro-trame et de la synchro-bit est effectuee par correlation sur Ia phase 
differentielle. 

Le signal complexe normalise retranspose en bande de base & la reception a pour expression : 
S(t) = exp{j[2*Af 0 t + do + *(t)} 
30 dans laquelle Afo represente I'ecart entre ia frequence d'emission f e et la frequence de reception f r auquel 
s'ajoute la frequence d'effet Doppler f d lorsque le recepteur se trouve en mouvement par rapport a 
I'emetteur : 

Afo = f e - fr + fd- 

Grace a la sequence preiiminaire, on dispose d'un signal de reference connu sur une duree NT, soit : 
35 R(t) = exp{j<fr(t)} avec te[0,NT] 

On definit alors un signal s' egal au produit du signal S par son conjugue* retarde* de deux temps bit : 
S(t) = S(t).S-(t-2T) = exp{j[4irAfoT + <p{t) - <f>(t-2T)]} 

En posant A<|>(t) = <*>(t)-<fr(t-2T) on fait apparattre la phase differentielle entre deux temps bit, 
d'ou S'(t) - exp{j[4*rAf 0 T + Af(t)]} 
40 Dans cette expression de S (t), le terme de phase initiate a disparu et I'ecart de frequence se traduit par 
un dephasage constant. 

On definit de meme a partir du signal de reference R un autre signal R tel que : 
R(t) = R(t).R*(t-2T) - exp{jA<>(t)> 

La fonction de correlation des deux signaux complexes S' et R' s'ecrit : 

45 

NT 

C(t) = J S' (tJ.R'Mt-rJdt = 
T 

NT 

50 exp{j4MfoT} T J exp{j[Acp(t)-A«)(t«T)]}dt 



En prenant le module au carre de C(t), le terme exp{j47rAf 0 T} disparaTt: 



NT 



|C(T)|2 = | / exp{jU<p(t)-Atp(t-T)]}dt!2 
T 
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La recherche du maximum de la fonction |C(t)F permet alors de determiner le debut du paquet. parce 
que |C(t)F est maximum pour r = 0. 

L'intgret d'effectuer la correlation sur la phase differentielle est que le pic de correlation est independant 
de la phase k Porigine et peu dependant de Tecart en frequence tant que Af 0 T«1. c'est-&-dire tant qu'on 
5 se situe dans la bande de transmission du filtre plac6 en sortie de I'gmetteur. 

Cependant le niveau du maximum du pic de correlation est plus sensible au bruit (ce qui revient a avoir 
une degradation de 3 dB du rapport signal k baiH). 

Deux aspects interviennent dans le choix de la sequence preliminaire : sa longueur (N = nombre de 
bits) et la configuration des bits. 
io Plus la sequence sera longue, meilleures seront les Probability de Fausse Alanme (PFA) et de Non- 
DStection (PUD). 

La configuration binaire de la sequence a une influence sur la precision du calage temporel. Le choix 
n'est pas trfcs facile, mais on peut malgre tout la choisir en respectant les contraintes suivantes : 
• suite non p^riodique (sinon il se forme plusieurs pics de correlation), 
75 • suite non constante (sinon il apparattun etalement temporel important), 

- suite qui ne comporte pas trop de valeurs alternees (sinon cela entraTne des variations de phase trop 
faibles). 

On peut appliquer la methode de correlation exposee ci-dessus au calcul de ia synchro-trame en 
effectuant cette correlation a raison de 2 <§chantillons par temps bit. En designant la valeur S(i T/2) de S(t) 
20 pour t = i T/2 et i entier par S(i) on a : 
S(i) = exp{j[27rAfoi I + Oa + z )]} 
Rifi) = exp{j>(i2)} 
On calcule : m 

S'(i) - Sfi).S ( i-4) = exp{j[4irAf 0 T + A$(i \ )]} 
25 R(i) = R1CO.R1* (i-4) = exp{jA*(i \ )} 
En posant : 
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Ct(j) = | J S'(i).R'i(i-j)jf 

la recherche du maximum de CiQ) permet de detecter le debut du paquet Lorsque ce maximum est 
detecte, la synchro-trame est acquise et I'instant d'emission du paquet est connu a t T/4. 
Cette precision n'est pas suffisante pour determiner la synchro-bit. 

Pour affiner I'estimation, il faut effectuer une seconde correlation avec un signal de reference decale de 
T/4. On dSfinit comme pour la premiere correlation : 
R 2 (i) = exp[j*<i T/2 + T/4)] 
R 2(i) = expOA*(i T/2 + T/4)] 

c 2 (j) = | 2 |~V(i)R'J(i-j)| 2 

^ Les deux correlations Ci(j) et C^O) vont presenter des pics pour les indices Ji et J 2 respectivement. 
Si Ci(Ji)£C2(J2) on prendra synchro-bit = J1 T/2 
Si C2(J2)>Ci(Ji) on prendra synchro-bit = J2 T/2 + T/4 

Cette double correlation permet d'avoir une precision de ± T/8 sur la synchro-bit. 
Compte tenu de la lenteur de la variation de la phase pour la modulation GMSK, cette precision est 
so suffisamment bonne pour connaTtre I'instant d'Schantilionnage. 

On a d6ja precise que dans certains cas 0C1 le rapport signal k bruit est faible (Eb/N 0 <6 dB), le choix 
de la synchro-bit pouvait §tre erronS. Cela resulte de I'appreciation incertaine du niveau des pics de 
correlation (par exemple si Ton choisit Jz T/2 + T/4 au lieu de J1 T/2). 

Pour faire rultime decision on va dSfinir SYNP la vaJeur de synchro-bit dite "principaie" et SYNS la 
55 valeur de synchro-bit dite "secondare". 
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Si Ci(Ji)SC 2 (J2) 


alors 


SYNP: 


Ji T/2 






SYNS : 


J 2 T/2 + T/4 


Si C 2 (J2)>Ci(Ji) 


alors 


SYNP: 


J 2 T/2 + T/4 






SYNS: 


Ji T/2 
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Le filtrage adapte sera effectue avec la valeur SYNP. 

La modification eventuelle sera effectuee dans I'algorithme d'estimation de I'ecart de frequence et de la 
phase initiate a partir d'un critere d'erreur qui sera defini ulterieurement. 

L'erreur sur la synchro-bit a peu d'effet sur le filtrage adapte ; par contre cela degrade enormement 
I'estimation de Afo et So. 

Cette methode de double correlation sur la phase differentielle est judicieuse car elle permet de 
determiner le debut de la trame et d'effectuer une premiere estimation de la synchro-bit. Ce dernier 
parametre sera confirme ou ajuste au cours de I'estimation de la phase de la porteuse. 



15 
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Filtrage adapte 

On peut montrer (cf. P.A.LAURENT : "Interpretation des modulations d'indice demi-entier. Extension a 
des indices voisins et applications", 9eme Coiloque GRETSI, Nice, Mai 1983, pages 503 a 509) que toutes 
les modulations numeriques de la forme S(t,B) = exp.j.*(t,B) peuvent se representor sous la forme d'une 
modulation d'amplitude selon I'expression suivante : 



25 



S(t,B) = . I 



j^F^t-iT) 



n 

avec C.= IT B. 

1 j,3S-o* i 
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dans laquelle F p (t) est la fonction principale. 

La decomposition de la modulation GMSK en une modulation d'amplitude est particulierement interes- 
sante, car elle permet de determiner facilement le filtre adapte. 

Ce dernier a une reponse impulsionnelle egale a F p (t - synchro-bit). 

Le terme synchro-bit tient compte de la position du signal regu par rapport a I'horloge d'echantillonna- 



ge. 



Le filtre adapte est realise sous la forme d'un filtre a Reponse Impulsionnelle Finie a 1 1 coefficients. 
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Estimation de la phase initiate et de I'ecart de frequence 

La methode envisagee est basee sur Sexploitation de la sequence preliminaire. 

En sortie du filtre adapte on dispose d'un signal regu dont la variation en fonction du temps est 
representee sur la figure 5a. Apres normalisation, ce signal a pour expression : 
Z(t) = exp{j[2-7rAfot + Bo + £(t)} 

La sequence preliminaire etant connue, on peut facilement calculer revolution du signal sur un intervaile 
[0,NT], N etant le nombre de bits de la sequence preliminaire. 

On connatt alors le signal de reference dont la figure 5b montre la variation en fonction du temps et 
dont I'expression normalisee peut s'ecrire : 
2o(t) = exp{j*(t)} pour te[0,NT] 

En effectuant le produit du signal regu Z(t) par le conjugue du signal de reference Z 0 (t) on eiimine le 
terme $(t) dO a la modulation (figure 5c). 
Z(t).Z 0 ( t) ■ exp{j[2^Af 0 t + Bo]}. 

L'etape suivante consiste a transformer le signal complexe obtenu en une variation lineaire traduisant 
revolution de la phase. 

0(t) =Arctg[Z(t).Z o (t)] - [2->rAfot + 0o]moduio Ztt 

Pour cela, it faut derouler la phase en eliminant les sauts de phase de 2*r. 

La figure 5d represente cette variation dont ('equation est : 
y = Aojo-x + Qo avec A«o - 2^ Ah . ^ 

Par une methode de regression lineaire, on peut calculer les parametres estimes A« 0 et ?o. Ce calcul 
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est systdmatique et done simple & mettre en osuvre. 

A partir des param&tres estimes, on peut maintenant comp8nser le signal en effectuant une multiplica- 
tion complexe : 



Signal regu : 
Signal compens6 



2(t) = exp{jtA«ot + $o + <Mt)]> 

S(t) = Z(t)expH(ASot+?o)] 

S(t) = exp{j[(Aa)o-A2o)t + eo-So + ^(t)]} 



10 Si l'estimation est correcte, 

aZo = = Bo et S(t) = exp{j4>(t)} 

Le signal obtenu n'est plus affecte par un ecart de frequence ni par la phase k I'origine. 
L'estimation de A«o et Bo est sensible k trois parametres : le bruit, !e calage tempore! et la longueur de 
la sequence preliminaire. 

'5 Lorsque le bruit augmente (Eb/N 0 <6 dB), cela peut entraTnerdes variations soudaines de phase qui se 
traduisent par des sauts de 2 n sur la phase deroulee. Ce probfcme a ete elimine en utilisant une technique 
de detection et de correction de saut de phase de 2 

La sensibility au calage tempore! est liee a revaluation de la synchro-bit. Si ce parametre est mal 
estime, la modulation n'est pas parfaitement eliminee ; cela se tradurt par une phase deroulee affectee d'un 
20 r^sidu de modulation. L'estimation de A** et Bo s'en trouve degradee. 

Pour avoir une estimation de frequence suffisamment bonne (erreur inferieure k 10 Hz) qui n'entratne 
pas d'erreur de decision, il est necessaire d'utiliser une sequence preliminaire de longueur superieure ou 
ggale a 64 bits. 

Sur 128 bits, cela entratneun rendement de transmission maximum de 50 %. 
25 Un tei rendement est tout k fait incompatible avec une transmission par paquets. 

Si Ton adopte au depart une sequence preliminaire de longueur plus faible N = 16 ou 32 bits, la 
methode exposSe ci-dessus permet d'obtenir des echantillons compenses. 

Cependant, la precision de l'estimation n'est pas suffisante pour corriger parfaitement la phase lorsque 
le niveau de bruit est important. 
30 une erreur de 20 Hz entre le debut du message et la fin se traduit par une rotation de phase de 58 , 
ce qui entratnedes erreurs de decision sur la fin du paquet. 

Ce sont done les bits proches de la fin du paquet qui sont les plus affectes. 

L'idee de I'invention consiste a decider un certain nombre de bits, par exemple les 16 bite suivant la 
sequence preliminaire et de refaire le processus d'estimation en considerant une nouvelle sequence de 
35 reference correspondant aux N bits de sequence preliminaire plus 16 nouveaux bits decides. 

En quatre nouvelles passes, on arrive a obtenir une precision de quelques Hertz pour Eb/N 0 = 6 dB. 
Cela fait a la fin une estimation sur N + 64 bits. 

La longueur de la sequence preliminaire a pu etre ainsi notablement reduite par cette methode 
d'estimation en plusieurs passes qui exploite les decisions intermediates sur des blocs de bits, qui resiste 
tres bien ku bruit et dont la convergence est rapide. 

Toutefois. ainsi qu'il a de\k &t6 evoque. l'estimation de A«o et Bo est sensible a la synchro-bit 

Une mauvaise estimation de synchro-bit va se traduire par une difference importante entre les points 
correspondant a la phase deroulee et la droite de regression : 
Soit 

45 

N | . - I2 

e = l y i~ (Au, °' x i + 9o) l 



Dans ce cas, € va augmenter de plus en plus. 

Trfcs rapidement (k la premiere ou a la deuxieme passe) e va depasser une vaieur de seuil et va 
ordonner un changement de synchro-bit 

On prendra alors pour vaieur de synchro-bit la vaieur secondaire SYNS. 

Tout le processus de caJcul est alors reinitialise pour refaire le filtrage adapte et l'estimation de A<*j et 
0o avec la nouvelle vaieur de synchro-bit 



8 



BNSDOCID: <EP. 



0349064A1J_> 



EP 0 349 064 A1 




It 



Decision 

Apr§s compensation, la decision est effectuee sur I'expression du signal mis sous la forme d'une 
modulation d'amplitude faisant intervenir la fonction principale F p (t). 
s Pour calculer les bits transmis il suffrt de faire enfin un d§codage differentiel. 

La figure 6 foumit sous forme d'organigramme tout !e processus de calcul. 

Le deroulement d'estimation des parametres Ao* et Bo & partir du DEBUT de programme (case 22) 
peut parattre lourd a mettre en oeuvre, mais ii est systematique et reiativement simple. 

On effectue d'abord une estimation grossiSre de la synchro-trame (case 23) et de la synchro-bit (case 
io 24) par correlation, calcul des vaieurs SYNP et SYNS pour la synchro-bit, et le choix initial SYN = SYNP 
pour la valeur de la synchro-bit SYN, ensuite le filtrage adapts (case 25), puis la determination approchee 
de la phase de la porteuse par regression lineaire sur 16 ou 32 temps bit : suppression de la modulation 
(case 26), deroulement de la phase (case 27), estimation de A«o et $o et calcul de I'ecart e (case 28). On 
compare enfin recart € calcuie a une valeur de seuil S (case 29), et on decide s'H agit de la dernfere passe 
75 (case 30). 

La suite du processus de demodulation peut alors etre decrite comme un systeme de deux boudes 
numeriques imbriquees : 

- une premiere boucle pour Pestimation de Ac*»o et So se refermant par la liaison 31, 

- une seconde boucle pour ('estimation de synchro-bit se refermant par la liaison 32. 

20 A chaque passage dans la premiere boucle num§rique, on red£cide dans la case 31 les N bits qui 
suivent la sequence preliminaire afin d'affiner au fur et & mesure I'estimation de A« 0 et So. Ces decisions 
sont dites intermediates. 

A chaque passage dans la seconde boucle num^rique, on calcuie dans la case 29 I'ecart e que Ton 
compare & la valeur de seuil. Ce critere va permettre de valider ou de reajuster la synchro-bit. Dans ce 
25 dernier cas, on prendra pour la valeur de la synchro-bit SYN la valeur secondare SYNS (case 32). 
On distingue deux cas possibles dans le traitement : 

- pas de remise en cause de la synchro-bit. La determination de Acjo et do est alors effectuee en quelques 
passes par une convergence rapide de la premiere boucle. 

- remise en cause de la synchro-bit. Dans ce cas, on recommence tout le processus de filtrage et de 
30 demodulation. Le temps de traitement par la voie de la seconde boucle devient alors plus long. 

La suite de I'organigramme a partir de la derniere passe (case 30) comporte la mise en oeuvre de la 
decision finale (case 33) et de la FIN du programme (case 34). 

A titre indicatif, on a represents sur la figure 7 les courbes A et B de taux d'erreurs (BER) theoriques 
relevees dans la litterature pour les modulations GMSK et MSK. 
35 Pour la modulation GMSK, les courbes 1 et 2 correspondent h un demodulateur coherent avec 

recuperation de porteuse par une boucle d'asservissement analogique, de bandes passantes B L = 460 Hz 
et 920 Hz respectivement 

Ce type de demodulation ne pourratt done pas fonctionner, ni en EVF ni en AMRT. De plus ces 
courbes sont obtenues sans ecart de frequence. 
40 La courbe 3 correspond au resultat obtenu avec un demodulateur differentiel analogique. 

A 10~ 2 de taux d'erreurs, la degradation par rapport a latheorie est tres importante (7 dB environ). 
Pour la modulation GMSK, les courbes de taux d'erreurs obtenus avec le precede de demodulation de 
invention sont portees sur les figures 8 et 9 pour des sequences preiiminaires de 32 et 16 bits 
respectivement et avec des ecarts de frequences Af 0 = 800 Hz (courbes 1) et Afo = 1600 Hz (courbes 2). 
45 Pour Afo - 800 Hz les resultats sont assez bons. 

A un taux d'erreur de 10~ 2 , on a une degradation de 1,2 dB pour N = 32 bits et de 1,4 dB pour N = 
1 6 bits par rapport aux taux d'erreur theoriques dont les courbes A et B (deja representees sur la figure 7) 
sont egalement reportees sur iesdites figures. 

Les resultats sont peu sensibles h I'ecart de frequence tant que Af 0 <l000 Hz. Au-dela t les resultats 
so sont un peu degrades. 

Le precede de Tinvention a permis de mettre en oeuvre un algorithme de demodulation coherente de la 
modulation numerique de type GMSK. 

Les resultats de simulation montrent que cette methode resiste bien au bruit et a un ecart de frequence 
m§me important entre Temetteur et le recepteur. 
55 Ce procede est done tout a fait compatible avec un fonctionnement en AMRT ou en EVF et il peut 

s'appliquer a n'importe quelle modulation presentant de I'interf^rence inter-symboie. 
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Revendications 



1. Procedg de demodulation coherente par traitement numerique d'un signal moduli numeriquement en 
phase continue et & enveloppe constants, le terms moduli de ladite phase Stant egal au produit de 
convolution de I'impulsion de phase etal(§e sur piusieurs temps bits par rinformation binaire transmise par 
paquets, le signal re$u transpose en bande de base sur deux voies en quadrature etant convert! en 
numerique et transmis h un processeur de signal qui effectue le traitement du processus de demodulation, 
caracteris<§ en ce que chacun desdits paquets comporte une sequence preliminaire par laquelie on dispose 
d'un signal de reference connu sur N bits permettant de detecter approximativement la synchro-trame et la 
synchro-bit par correlation avec la phase difterentielle et d'amorcer restimation des param&tres de phase 
initials flo et d'Scart residuel en frequence Afo, I'affinement progressif des valeurs approves Stant obtenu 
au moyen de deux boucles num<§riques imbriquees : une boucle lente initialise au-defc d'une valeur de 
seuil pour la detection de la synchro-bit, une boucle rapide effectuant des decisions intermtkilaires sur des 
blocs additionnels de bits pour restimation de do et Afo. 

2. ProcedS selon la revendication 1 , caracterise en ce que lesdites detections de synchro-trame et de 
synchro-bit sont obtenues par une premiere correlation sur la phase differentielle permettant de connaTtre 
rinstant d'emission desdits paquets k ± T/4, T etant la dun§e d'un bit, puis par une deuxieme correlation sur 
la phase differentielle effectu<§e avec un signal de Irenes d6cal<§ de T/4, les fonctions de correlation 
correspondantes prtsentant chacune un pic indgpendant de la phase initiate et tres peu dependant de 
Pdcart residuel en frequence, le pic de niveau supSrleur et le pic de niveau inferieur definissant respective- 
ment une synchro-bit principal SYNP et une synchro-bit secondaire SYNS, la precision de ± T/8 ainsi 
obtenue sur la synchro-bit 6tant suffisante pour connaTtre I'instant d'echantillonnage, ladite detection de 
synchro-bit etant suivie d'un filtrage adapts effectue avec la valeur SYNP par un filtre a reponse 
impulsionnelie finie de type Gaussian afln de limiter ta bande de bruit 

3. Proc§d6 selon les revendications 1 et 2. caracterise en ce que lesdites estimations de la phase 
initiate do et de I'Scart residuel en frequence Afo k la suite dudit filtrage comportent les etapes suivantes : 

- Elimination du terme de modulation en effectuant le produit du signal regu par le conjugue du signal de 
r£f£rencs 

- dSroulement de la phase en Sliminant les sauts de phase de 2* pour obtenir une variation liralaire ayant 
pour Equation y = Auox + So avec A«o = 2?rAfo, 

- calcur des paramfctres estimes A^o + ?o par une methode de regression iineaire et de I'ecart 



c = (Aw 0 .X.+ 9 0 ) 

1=0 11 1 

entre les points correspondant h ladite phase dEroulee et ladite droite de regression. 

4. ProcEde selon les revendications 1 a 3, caracterise en ce que ledit ecart e etant inferieur k ladite 
valeur de seuil, ladite estimation de Afo et So est affinee suivant une boucle rapide en piusieurs passes 
exploitant les decisions intermediates sur les N bits de la sequence preliminaire auxquels on ajoute a 
chaque passe un certain nombre de bits decidfes. 

5. Precede selon les revendications 1 a 3, caracterise en ce que ledit ecart € depassant ladite valeur de 
seuil du fait d'une evaluation d§fectueuse de synchro-bit ledit processus de calcul est reinitialise suivant 
une boucle lente pour refaire le filtrage adapts et restimation de Afo et 6o a partir de I'autre valeur de 
synchro-bit egale a ladite valeur secondaire SYNS. 

6. Procede selon la revendication 4, caracterise en ce qu'a la suite de la demiere passe, une 
compensation ne laisse subsister que la composante de phase du signal qui n'est plus affectee par I'ecart 
residuel en frequence ni par la phase h I'origine, la decision finale etant ensuite effectuee et un decodage 
differentiel foumissant enfin la suite d'informations binaires transmises. 

7. Procede selon I'ensemble des revendications 1 a 6, caracterise en ce qu'il s'applique aux 
demodulations coherentes de signaux modules selon des modulations de types GMSK, 2SRC, TFM, 
GTFM,... dont la loi devolution de la phase suit une variation progressive. 
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